home *** CD-ROM | disk | FTP | other *** search
/ WINMX Assorted Textfiles / Ebooks.tar / Text - Mathematics - Numerical Mathematics and Computing (F).zip / bisect1.f < prev    next >
Text File  |  2002-06-11  |  946b  |  46 lines
  1. C*****************************************************************************
  2. C
  3. C PAGE 77-78: NUMERICAL MATHEMATICS AND COMPUTING, CHENEY/KINCAID, 1985
  4. C
  5. C FILE: BISECT1.FOR
  6. C
  7. C FIRST VERSION OF BISECTION METHOD (BISECT,F,G)
  8. C
  9.       EXTERNAL F,G
  10.       DATA AF/0.0/, BF/1.0/, N/50/
  11.       DATA AG/0.5/, BG/2.0/ 
  12.       CALL BISECT(F,AF,BF,N)
  13.       CALL BISECT(G,AG,BG,N)
  14.       STOP
  15.       END 
  16.   
  17.       FUNCTION F(X) 
  18.       F = EXP(X) - 3.0*X
  19.       RETURN
  20.       END 
  21.   
  22.       FUNCTION G(T) 
  23.       G = T**3 - 2.0*SIN(T) 
  24.       RETURN
  25.       END 
  26.   
  27.       SUBROUTINE BISECT(F,A,B,N)
  28.       U=F(A)      
  29.       V=F(B)      
  30.       PRINT 6,A,U,B,V       
  31.       IF(U*V) 2,7,7 
  32.   2   DO 5 J=1,N  
  33.         C=(A+B)*0.5 
  34.         W=F(C)    
  35.         PRINT 6,C,W 
  36.         IF(W*U) 3,7,4       
  37.   3     B=C       
  38.         V=W       
  39.         GO TO 5   
  40.   4     A=C       
  41.         U=W       
  42.   5   CONTINUE    
  43.   6   FORMAT(2X,2E22.14)    
  44.   7   RETURN      
  45.       END
  46.